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Abstract 

The influence of the convective structure of the solar photosphere on the shifts of 
spectral lines of iron was studied. Line profiles in the visible and infrared spectrum were 
synthesized with the use of 2-D time-dependent hydrodynamic solar model atmospheres. 
The dependence of line shifts on excitation potential, wavelength, and line strength was 
analyzed, along with the depression contribution functions. The line shifts were found to 
depend on the location of the line formation region in convective cells and the difference 
between the line depression contributions from granules and intergranular lanes. In visible 
spectrum the weak and moderate lines are formed deep in the photosphere. Their effective 
line formation region is located in the central parts of granules, which make the major 
contribution to the absorption of spatially unresolved lines. The cores of strong lines 
are formed in upper photospheric layers where is formed reversed granulation due to 
convection reversal and physical conditions change drastically there. As a consequence 
the depression contributions in the strong line from intergranular lanes with downflows 
substantially increase. This accounts for smaller blue shifts of strong lines. In infrared 
spectrum the observed decrease in the blue line shifts is explained by the fact that their 
effective line formation regions lie higher in the photosphere and extend much further into 
the reversed granulation region due to the opacity rise with the increase of line wavelength. 
Additionally the effective line formation depths of the synthesized visible and infrared Fe I 
lines and their dependence on line parameters is discussed. 

1 Introduction 

If the motion of the Sun-Earth system is taken into account as well as the correction for the 
gravitational red shift (636 m/s) caused by the difference between the gravitational potential 
of the solar atmosphere and the Earth is introduced, the most absorption lines observed at the 
center of the solar disk are shifted towards shorter wavelengths with respect to laboratory lines. 
The observed blueshifts of the solar iron lines are 300-400 m/s on the average [H El HI [5l El 
[271 [29l [33] . They were measured to an accuracy of 100 m/s in the range from to 1 km/s. It 
is found that the blueshifts decrease with growing height of line formation as suggested by its 
dependence on excitation potential. The largest blueshifts are observed in weak lines with high 
excitation potentials, while the smallest blueshifts are in strong lines. The redshifts were found 
to be small in some strong lines in the visible spectrum, but they can be as large as 1 km/s in the 
ultraviolet lines in the range 195-200 nm [33]. Some "reddening" of blueshifts was also noted in 
active photospheric regions [3l H] . The blueshifts observed in weak and moderate lines can be 
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explained by a decrease with height in the velocity of the convective motions along the line of 
sight towards the observer ^ [B] . The formation region of strong photospheric lines extends to 
the layers of reversed granulation, a pattern comparable to the granulation but with reversed 
brightness modulation. The matter at the center of convection cells becomes cooler than in 
descending intergranular lanes, but it still moves upwards. An inverse correlation between 
brightness and shifts was found in the cores of very strong lines in the visible and infrared 
spectrum [23]. The effects of convection overshoot in the solar photosphere [28] is believed to 
be the cause of the observed line redshifts. 

The objective of this study is the effect of reversed granulation and the role of granules and 
intergranular lanes in the formation of spatially unresolved spectral lines and the cause of the 
"reddening" of blueshifts of strong lines. Such a study is possible today due to the remarkable 
progress in the simulation of solar granulation and its evolution [21 [TJ |28l |3T1 |3l] and in the 
development of the idea of depression contribution functions [201 1211 1211 [2S]- Based on 2-D 
hydrodynamic model atmospheres [TH], we synthesized Fe I lines and calculated their shifts. 
We used contribution functions to investigate the region of line formation in convection cells 
and to analyze the influence of physical conditions in the granular photosphere on the shifts of 
the synthesized lines. 

2 Hydrodynamic model of solar granulation 

The sequence of single-scale (one-cell) and multi-scale two-dimensional nonstationary hydro- 
dynamic (HD) models of solar granulation built by Gadun [H [19] were used more than once 
for the synthesis of photospheric lines [3 UHl IHI 1121 1131 HH HSl 113 HB] ■ Such sequence models 
proved to be capable of reproducing satisfactorily the manifestations of convective motions, the 
asymmetry and shifts of lines in particular. Comparison of the profiles synthesized with the 
use of 3-D and 2-D HD models [8j showed that the single-scale 2-D models are inferior to 3-D 
models in reproducing the observed line asymmetry because their inhomogeneity spectrum is 
narrower. Nevertheless, the use of 2-D models [TSj for the synthesis of spectral lines is justified 
when we want to obtain the overall picture of the processes in the superadiabatic regions in 
the photospheres of the Sun and solar-like stars. 

The 2-D HD model sequence we used in our study treats the radiative transfer on the 
basis of the monochromatic absorption coefficient depending on frequency, and the atomic and 
molecular line absorption is taken into account through the use of the opacity distribution 
function (ODF) method [22]. The thermal convection in the solar envelope was assumed to be 
quasi-stationary with a single convective flow in the simulation region, and this flow represented 
the only granulation scale. In building the HD models we assumed that convective flows were 
quasi-stationary in time, the horizontal size of convective cells did not change in the course 
of their evolution, and convective cylinders always extended downwards as deep as permitted 
by the simulation region (1100 km below the surface level in the models). These models were 
called the steady-stable (SS) models by their author. The dimensions of the calculated region 
were 1400 km in the horizontal direction and 1960 km in the vertical direction (the atmospheric 
layers extended to 800 km). The spatial step was 28 km. The total duration of simulation was 
50 min of the real solar time with a time step of 30 s between the models. Test calculations 
showed that such a duration was sufficient for a correct representation of five-minute oscillations 
in the spectral line calculations. This time-dependent sequence of 2-D models for the solar disk 
center includes 100 models (or 100 snapshots) consisting of 48 columns. Each column is located 
along the line of sight. A more detailed description of 2-D HD models of this type can be found 
in [91 [19]. 
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3 Spectral line synthesis 



To find the dependence of line sliifts on principal line parameters, we used certain artificial Fe I 
lines. The central line depths (d), lower excitation potentials (EP), and wavelengths of the 
artificial lines was specified in a wide range of actual Fe I line parameters. Then the product of 
abundance and oscillator strength Agf was chosen for every line in accordance with its central 
depth. The lines were synthesized in the 1.5-D LTE solution of the transfer equation [TB] with 
the use of the SS-model sequence - this means that the line profiles were calculated for each 
column in all snapshots as for a plane-parallel atmosphere. Then the obtained line profiles were 
averaged over all columns in each snapshot and over all snapshots. The result was a line profile 
averaged over space and time. 

The central wavelength of a synthetic profile A was determined by fitting a fourth-degree 
polynomial to the profile core. The accuracy of the procedure was ±50 m/s. The shift of 
the synthetic line with respect to the laboratory central wavelength Aq was calculated, with 
allowance made for the gravitational shift, in radial velocity units by the formula Vr = c(A — 
Ao)/Ao — 2.12 ■ 10~^c, c being the light velocity. Test calculations demonstrated that the shift 
of the averaged line profile synthesized with the full SS-model sequence (100 snapshots) differs 
insignificantly from the shift of the line profile calculated with a snapshot averaged over all 100 
snapshots. This allowed us to substantially cut the calculations by using the averaged SS model 
only. For line shift analysis we used the calculated line profiles for each column in the averaged 
SS model and the profile averaged over all columns. We also used the depression contribution 
functions calculated for each profile point and the effective formation heights. These latter are 
calculated as the weighted mean for the depression contribution function [201 EI] ■ 

4 Analysis of calculation result for lines in the visible 
spectrum 

The plots of line shifts of synthetic Fe I lines for various wavelengths, line strengths, and heights 
of formation (Fig. 1) suggest that these calculated shifts behave in the same fashion that the 
observed shifts [H El |6]. Weak lines with high EP at shorter wavelengths have the largest 
blueshifts that correspond to the largest negative radial velocities. Strong lines, whose central 
parts are formed high in the photosphere, feature redshifts. The typical strong narrow lines 
with d = 0.8, equivalent width W = 8 pm, and A = 600 nm show the largest redshifts (positive 
radial velocities). Our results also demonstrate a clear dependence of shifts on wavelength 
(Fig. la), although no such dependence was found in the observed line shifts [1]. It seems 
likely that the authors could not reveal this dependence because they ignored the differences in 
the central depths of the lines observed. In general, we confirm that the line blueshifts decrease 
with growing line wavelength, line strength, and effective height of line formation. 

To understand the causes of formation line shifts, we considered the snapshot of the sim- 
ulated convective cell (the averaged SS 2-D model) and the depression contribution functions 
[CF) of synthetic Fe I lines with different line parameters on Fig. 2. The isotherms (dashed 
lines), which demonstrate the temperature distribution in the convective cell, show that the 
temperature inversion in the photospheric layers begins above 200 km. This is the region of 
the so-called reversed granulation caused by convective reversal. Here the isotherms, which are 
convex over the granule center, become concave. The velocity field marked by arrows illustrates 
the motions in the convective cell: vertical upfiows in the central part and downfiows at the 
cell periphery. The central region of convective cell will be call granule, and the regions with 
downfiows will be call intergranular lanes, as adopted for the observed granulation structure in 
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Figure 1: Dependence of convective shifts of synthesized Fe I hnes on wavelength (a), central 
depth (b), equivalent width (c), and height of formation (d). Solid line: EP = 2.5 eV, dotted 
line: EP = 4.25 eV, dashed line: EP = 6.4 eV. Size of the squares grows with line strength 
(a), d = 0.1, 0.3, 0.5, 0.6, 0.8, and excitation potential (b,c,d), EP = 2.5, 4.25, 6.4 eV. 
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the photosphere. In addition, we also distinguish a region with predominantly horizontal flows 
of matter located between the granule and intergranular lanes. The depression contribution 
function contours calculated in each column for the core of synthesized line were plotted on 
the snapshot of the convective cell from the averaged SS 2-D model. The solid thick line is 
the contour corresponding 0.5CFmax- It outlines the region of effective depression contribu- 
tions. Inside this contour, there are the contour of the greatest contributions corresponding 
0.75CFmax, and outside it two contours corresponding 0.25CFmax and O.lCFmax are plotted. 
Here CF^ax is the maximum of contribution functions calculated in each model column for 
line core. Outside the outermost contour the line absorption is insignificant. Fig 2 shows the 
contours for the weak line (a) and strong line (b); for moderate line with EP = 2.5 (c) and 
EP = 4.25 (d). Figs 2e,f,g,h show the CF contours for moderately strong line depending on 
the wavelength. It is interesting that lines with wavelengths about A 800 nm are only partially 
formed in the granulation inversion zone, while the lines with AA 1600-1800 nm are almost 
completely formed there. The depression contribution from intergranular lanes is greater for 
infrared lines. Figure 2h illustrates the appearance of new small effective contribution regions in 
intergranular lanes. In general, the depression contribution contours for different lines plotted 
on the background of a snapshot of a convection cell clearly demonstrated that intergranular 
lanes can play different part in the formation of absorption lines depending on line wavelength 
and central line depth. 

Apart from contribution functions, we also considered the individual line profiles calculated 
for each column in the simulated convective cell. We selected the most typical line profiles 
formed in intergranular lanes, in region with predominantly horizontal motions, and in the 
central parts of the granule. Figure 3 displays these line profiles as well as the profiles of the 
line formation heights. It should be noted that the line formation height profiles are quite 
intricate for inhomogeneous photosphere models as different from the homogeneous models. In 
the latter case these height profiles are always completely symmetric. The synthetic line profiles 
at the center of granule are slightly asymmetric and have always blueshifts (kind I). The profiles 
in the regions with predominantly horizontal motions have weak asymmetry and are virtually 
not shifted (kind II). They differ insignificantly from kind I profiles in strength and width. The 
profiles in intergranular lanes are strongly asymmetric and have large redshifts (kind III). Their 
cores are weaker and they are markedly widened as compared to the profiles of kinds I and II. 
Depending on the line parameters the shape and shift of the profile of I, II, III kinds slightly 
change, but its major features remain the same. It should be noted that the line profiles observed 
in the spectra obtained with very high spatial and spectral resolution for bright, moderately 
bright, and dark regions of photospheric granulation P2] are in complete agreement with our 
synthetic profiles. This confirms to a degree the correctness of our calculations with the use of 
2-D HD models. 

Comparison of the profile averaged over the whole simulation region with the profiles of 
kinds I, II, and III (Fig. 3a) shows that the largest contribution to the averaged profile core 
of a weak line {d = 0.1) comes from kind I and II profiles, i.e., from the whole granule and 
region with predominantly horizontal motions. The contribution of the kind III profiles to the 
average profile of the weak line is noticeable in the red wing of the line only — the wing is more 
extended, increasing the profile asymmetry in extended wings. The blue wing of the average 
profile of the weak line, as the line core, is formed in the main by kind I and II profiles. It 
is therefore concluded that the core shift and the shape of the blue wing of the weak line are 
controlled primarily by the ascending motions in the central parts of granules, while the shape 
of the red wing is also affected by the descending motions in intergranular lanes. From Fig. 3c 
we can estimate the effective heights of formation of the core and the wings in the averaged 
profile of a weak line [d = 0.1). This is the layers in the range of 140-160 km for the core, 80- 
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Figure 2: The snapshot of the convective cell from the averaged SS 2-D model and contours 
of depression contribution functions (solid line) for various synthesized Fe I lines: weak line 
and strong line (a,b); moderate line with low and high excitation potential (c,d); moderate line 
with different wavelengths (e,f,g,h). Thick solid line marks contour corresponding 0.5CFmax- 
Dashed line: isotherms (4000, 4000, 6000, 7000, 8000, 9000, 10000 K from the top down), sohd 
dashed line: surface level with the Rosseland optical depth equal to unity. Arrows indicate the 
velocity direction and their length is proportional to the velocity amplitude. 
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Figure 3: Three kinds of line profiles calculated in the simulated convection cell for a weak line 
(a) and a strong line (b) at A = 500 nm, EP = 4.25 eV. Dot-and-dash line: kind I profile from 
the central region, dashed line: kind II profile from the region with horizontal motions, dotted 
line: kind III profile from intergranular lanes; solid line: average profile. The profile of effective 
heights of line formation (c,d) are in response to changing line profiles (a,b). 

140 km the blue wing, and 50-140 km for the red wing. The effective heights of the moderate 
line {d = 0.5) is located 150 km higher then the above weak line. All other details relative to 
different kinds of profiles are similar. In the strong line {d = 0.8, Fig. 3b) the kind III profiles 
formed in intergranular lanes are much stronger, and the line asymmetry is stronger because of 
temperature drop in the upper granular layers. As a result, the shape of the unresolved profile 
of the strong line change, and the blueshift of the line core decreases. The strong line core is 
formed at heights of 280-350 km (Fig. 3d) in the reversed granulation region. 

The above analysis of the shape and shifts of synthesized line profiles allows us to conclude 
that the red wing of spatially unresolved line may be regarded as an indicator of the role which 
intergranular lanes play in the formation of absorption lines. The more extended red wind as 
compared to the blue wing, the larger is the depression contribution in the spectral line from 
interganular lanes. The shape of blue wings and shifts of weak and moderate lines in particular, 
depend completely on the physical conditions in the central parts of granules. 

5 Analysis of infrared lines 

Lists of numerous unblended infrared lines in three bands of the solar spectrum - J (1.00-1.34 
/im), H (1.5-1.8 fj,m), and K (1.9-2.5 /im) - were published recently in [301 [32]. The H band 
lines are most frequently used in various spectral investigations of stellar atmospheres. These 
lines have some advantages over other infrared lines. The continuum opacity is at a minimum 
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in the H band, and wc can see deeper photospheric layers in these hnes as compared to the 
hnes of J and K bands or the visible lines. For instance, the Fe I lines with very high excita- 
tion potentials and large Lande factors are very sensitive to the magnetic fields located at the 
photosphere base. It is also well known that the contrast of continuum intensity is smaller for 
infrared lines, and this allows us to suggest that they are not formed in granulation cells in the 
same way as the lines in the visible range. Some differences between the profiles of the infrared 
and visible lines with the same principal parameters are attributed to the strong wavelength 
dependence of absorption coefficient. The selective absorption grows proportionally to wave- 
length, while the continuous absorption depends to a great extent on the sort of absorbers. 
Recall that the continuum opacity in the infrared range attains its maximum at A 1000 nm, 
diminishes with increasing wavelength to a minimum at A 1600 nm which is approximately 
of the same magnitude as the minimum at A 400 nm, and then grows again. Generally, the 
total opacity is always greater in the infrared lines than in the lines of the visible spectrum due 
to the increase of the selective absorption coefficient. As a result, the regions of formation of 
any infrared lines are located in higher layers and are more extended into intergranular lanes. 
These effects increase with line strength and excitation potential. The excitation potentials 
for the infrared Fe I lines range from 2 to 6.4 eV, and for a line with d = 0.5, A = 1800 nm, 
and EP = 4.25 eV, for instance, the region of its core formation is located completely in the 
reversed granulation layers. In intergranular lanes the line formation region shifts higher (Fig. 
2h) with respect to the visible lines of the same hne strength (Fig. 2e). The depression contri- 
bution from intergranular lanes in the core of a visible line is about 25 percent of the maximum 
contribution, while for infrared lines it is 50 percent and more. 

The depression contributions of intergranular lanes to the infrared line affects the shape of 
average line profiles. Figure 4 shows line profiles with a central line depth of 0.5 and wavelengths 
of 400, 800, 1600, and 1800 nm as well as their formation height profiles. All infrared line 
profiles are nearly twice as wide as the profiles of similar visible lines. The widest and the 
most asymmetric profile is formed in intergranular lanes, where the temperature and velocity 
gradients are large, and as a result the asymmetry in the red wings of average profiles is greater 
than in the visible lines. The effective height profiles also suggest that the formation regions of 
infrared lines are more extended in height. Weak lines are formed in the range from —100 km 
to 300 km, while strong lines arc formed at heights from 100 km to 330 km. 

Fig. 5 shows the effective heights of formation of Fe I lines in the infrared and visible 
spectrum. It should be noted that the weak iron lines {d < 0.1, AA 1600-1800 nm, EP = 6.4 eV) 
are formed in the deepest photospheric layers. The blueshifts of these deepest-formed infrared 
lines are not the largest ones (see Fig la) due to the effect of additional absorption from the 
lanes. Therefore the blueshifts of infrared weakest lines are smaller than the blueshifts of the 
visible weakest lines with EP = 4.25 eV. 

Thus, the infrared iron lines are mainly formed in the reversed granulation zone. Therefore 
intergranular lanes play considerable importance role in the formation of the cores of infrared 
lines. This accounts for the observed decrease of blueshifts of infrared spectrum as compared 
with the visible spectrum. 

6 Conclusion 

We synthesized the Fe I lines in the framework of the 2-D HD models of solar granulation and 
determined the wavelength shifts of these lines and their dependence on the line parameters. 
The derived relations fit quite well the observation data. Additionally we calculated the depres- 
sion contribution functions to demonstrate that the line shifts are determined by the location 
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Figure 4: The same as in Fig. 3, but for moderate lines with d — 0.5, EP — 4.25 eV, and A400 
(a), A800 (b), A1600 (e), A1800 nm (f). 
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Figure 5: Effective heights of formation for artificial synthesized Fe I lines depending on 
wavelength (a) and central depth for the visible (b) and infrared (c) spectrum. Sohd hue: 
EP = 2.5 eV, dotted line: EP = 4.25 eV, dashed line: EP = 6.4 eV. Size of the squares grows 
with line strength d = 0.1, 0.3, 0.5, 0.6, 0.8 (a), and wavelength A = 300, 400, 500, 600, 800 nm 
(b) and A = 1000, 1200, 1800 nm (c). Crosses mark A = 1600 nm. Dot-and-dash line marks 
A = 400 nm, EP = 4.25 eV shown for comparison. 
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of effective line formation region in convection cells. Obtained results confirmed the convective 
character of the line shifts. 

Analysis of the line profiles and the corresponding contribution functions showed that the 
contribution from granules and intergranular lanes to the depression of spatially unresolved 
spectral lines is essentially different for weak and strong lines as well as for lines in the visible and 
infrared ranges. In the visible range weak and moderate lines originate deep in the photosphere, 
where the classical convection structure is not disrupted, i.e., hot masses move upwards and 
cool masses descend. Cores of strong iron lines are formed in the upper photosphere, in the 
convection reversal zone, where occurs a temperature inversion as well as granulation inversion. 
The temperature in the ascending central regions of a convection cell is lower than in descending 
regions. The temperature inversion increases the role of intergranular lanes in the formation of 
the depression of the spatially averaged line. As a result, the depression contribution from the 
regions with downfiows grows, and the line blueshifts diminish to zero or even become red. In 
the infrared spectrum the effect of reserved granulation on the line formation is stronger. Due 
to growth of hue opacity with wavelength the level of effective hue formation is shifted higher 
into the photosphere where are manifested the reserved granulation zone. Therefore the line 
blueshifts in the infrared spectrum are smaller than in the visible spectrum. 

Convective motions cause the asymmetry of spectral line profiles. The blue wings and cores 
of spectral lines are mainly formed in the central parts of granules, where upfiows prevail. The 
red wings are mainly formed in intergranular lanes with downfiows. Greater contributions from 
intergranular lanes in the depression of strong lines results in the strengthening of red wings 
of absorption lines and consequently a greater asymmetry, as compared to weak and moderate 
lines. Therefore the red wing of a spectral line can be regarded as an indicator of the additional 
absorption in intergranular lanes. 

We have found that the weak infrared Fe I lines with A ^ 1600 nm and very high excitation 
potentials (about 6.4 eV) are formed at the deepest photospheric layers. Such lines may be of 
prime importance in diagnostics of stellar atmospheres. 
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